Please note this document is a working draft.
Cities have used urban green spaces as a tool to enhance the health of residents, increase biodiversity, and mitigate the effects of climate change. Research has shown that urban green space can alleviate urban-heat-island effects, improve air quality, and encourage outdoor socialization and physical activity.1 However, urban green spaces are more frequently supported in higher-income, predominantly white neighborhoods, exacerbating inequity in health.2 Several U.S. cities, including Denver, have proposed goals to equitably increase green space in their cities.
Often, urban planners and residents use turf and exotic plants to expand green space, which are are not ideal for pollinator health and demand high amounts of water, which may not be sustainable in the mountain west, a region expected to experience more severe drought due to climate change. While a growing body of research suggests exposure to urban green space can improve human health and prevent premature mortality,1,3 little is known about the potential impact of native vegetation on human health in cities in the U.S. mountain west, a region with a growing population vulnerable to the impacts of climate change.
We aim to estimate the number of premature deaths that would be prevented by the implementation of various policy scenarios adding green space, including native vegetation, to the City of Denver. These policy scenarios were informed by conversations with local stakeholders in advocacy, research, and governance.
We measured existing green space in the City of Denver by the normalized difference vegetation index (NDVI) measured by the Landsat-8 satellite at a spatial resolution of 30 square meters. To define native-plants greenness, we measured the NDVI of several places in Denver (places described in more detail here) comprised of various levels of native and adapted plants during cloud-free summer (May, June, July, August) days over five years (2016-2021). Over these days in these areas, the mean NDVI was as low as 0.35 at the Denver Botanic Gardens Green Roof and as high as 0.67 in City Park, which is 0% native and consists of more watered turf. A plot in Denver Botanic Gardens that is 100% native plants had a mean NDVI of 0.54 over these days.
| Place type | Place name | NDVI, mean | NDVI, 25th-ile | NDVI, 50th-ile | NDVI, 75th-ile |
|---|---|---|---|---|---|
| high diversity | Chatfield Meadow Restoration | 0.48 | 0.47 | 0.50 | 0.50 |
| high diversity | Chatfield Prairie Garden | 0.37 | 0.35 | 0.37 | 0.37 |
| high diversity | City Park Greenhouses | 0.31 | 0.27 | 0.28 | 0.28 |
| high diversity | Denver Botanic Gardens Green Roof | 0.34 | 0.34 | 0.36 | 0.36 |
| high diversity | Kenderick Lake Xeriscape Garden | 0.40 | 0.36 | 0.37 | 0.37 |
| native spectrum | Denver Botanic Gardens, 100% Native | 0.54 | 0.52 | 0.55 | 0.55 |
| native spectrum | Plains Conservation Center, 100% Native | 0.16 | 0.23 | 0.26 | 0.26 |
| native spectrum | Green Mountain Park, 85% Native | 0.42 | 0.36 | 0.40 | 0.40 |
| native spectrum | Hogback along C-470, 75% Native | 0.36 | 0.32 | 0.35 | 0.35 |
| native spectrum | Suburban Open Space 1, Chatfield H.S., 50% Native | 0.40 | 0.37 | 0.41 | 0.41 |
| native spectrum | Suburban Open Space 2, Columbine Hills Church, 30% Native | 0.36 | 0.26 | 0.31 | 0.31 |
| native spectrum | Denver Botanic Gardens at Chatfield, 10% Native | 0.53 | 0.40 | 0.54 | 0.54 |
| native spectrum | City Park, 0% Native | 0.67 | 0.67 | 0.70 | 0.70 |
To acknowledge the lack of precision in these NDVI estimates and the range of NDVI levels that may correspond to a native-plants intervention, we defined NDVI using the mean of the Denver Botanic Gardens as as the primary analysis and that of the Denver Botanic Gardens Green Roof as an alternative value.
We conducted ten interviews with local stakeholders who are actively engaged in advocacy, policy, and governance related to native plants and green space in Greater Denver. Following these interviews, we developed four policy scenarios informed by the collective set of visions and priorities articulated by the stakeholders. The scenarios range from short-term goals to more ambitious, perhaps less immediately realistic visions. The scenarios fall into five categories.
Homogenous native-plants greening: Under the first set of scenarios, we consider adding native plants homogenously across all block groups, without respect to the specific location. Specifically, we consider two scenarios:
Set 30% of the block group’s area to be as green as native plants, corresponding with a goal set by some scientists and advocates to protect 30% of lands and oceans by 20304;
Set 20% of the block group’s area to be as green as native plants, a level that, according to some of our interviewees is sufficient to support pollinator health.
Equity-driven native plants greening: Under the second set of scenarios, we propose adding native plants to the same amount of area as above (30% or 20% of Denver’s area) but rather than applying the intervention homogeneously across all black groups, we prioritize application of the native plants intervention to those areas of greatest need from a health-equity perspective. Both the State of Colorado and the City of Denver have created equity indices to guide their work. The Colorado Department of Public Health & Environment (CDPHE) describes their equity index here5 and makes the index available for download here.6 The City of Denver Department of Public Health & Environment (DDPHE) describes the neighborhood-level equity index here7 and makes the data available for download here.8 We use both definitions to consider a total of four scenarios:
Set 30% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the CDPHE definition of health equity.
Set 20% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the CDPHE definition of health equity.
Set 30% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the DDPHE definition of health equity.
Set 20% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the DDPHE definition of health equity.
Create native buffers around riparian areas (rivers, streams, lakes, and ponds) of the following sizes:
200 feet (ideal for ecosystem health but possibly not realistic)
100 feet (compromise)
50 feet (most realistic; less good for ecosystem health)
Initiatives related to green infrastructure and stormwater management. (Please refer to corresponding section for details.)
Replace a portion of parking-lot surface with native plants:
This section elaborates on the scenarios and visualizes them with maps.
Under the first scenario, which we view as the most ambitious, we consider the expected impact on mortality by applying a green-space intervention to all census block groups, without specifying where.
We first measured the mean NDVI of each census block group on July 4, 2021. The weighted block-group-level mean is presented here, weighted by the proportion each 30 square-meter pixel covers by the census block group. For example, if half of a pixel overlaps the block group, it receives a weight of 0.5 in the weighted average. We also removed bodies of water before measuring NDVI.
Figure 2.1: Scenario 1
For each definition of native-plants NDVI, We excluded census blocks whose baseline NDVI was above the threshold. This map visualizes the NDVI cutoffs both for lower of the two values, that of the Denver Botanic Gardens Green Roof, and the higher of the two, that of Denver Botanic Gardens Native Plot.
Under this set of scenarios, we assumed we had a “budget” of greening of either 30% or 20% of Denver’s area, and we applied that budget to those with most need from a health-equity standpoint based on the CDPHE and DDPHE definitions, described above. As the below maps show, the two health-equity definitions from the State and from the City overlap considerably.
We measured NDVI in a 200-foot buffer, a 100-foot buffer, and a 50-foot buffer around all bodies of water in Denver. We downloaded bodies of water from OpenStreetMap (code here). We defined residential exposure to these riparian areas as those individuals living within a 500-meter buffer, following research on green-space and health that have defined green-space exposure based on residential proximity.3 We estimated the number of people in this buffer by multiplying the population density of the census block group by the intersecting area.
The below map depicts mean NDVI in the portions of census block groups that intersect a 200-foot buffer as well as those pieces that intersect the part of the 500 m buffer that would not be intervened upon, i.e., the part between the 200-foot buffer and the edge of the 500 m buffer.
We spoke with representatives at the local Office of Green Infrastructure who described three categories of initiatives–some planned, some aspirational–that could add native or adapted plants to the city.
We were provided a list of planned projects throughout Denver. Per our conversations, we assume that about 75% of the project’s footprint would consist of native or vegetation. Like for riparian areas, we defined exposure to the projects as those individuals living within a 500 m buffer of the projects. A map of baseline NDVI of the projects themselves and of a 500 m buffer around the projects is below.
Not yet complete.
The approach for estimating the health impact of the redevelopment controls differs from scenarios 2 and 3.1. because we will not know where the re-developments will occur. We thus simulated possible locations. Per our conversations, we anticipate the following number of parcels will be subject to these rules per year:
We gathered data on parcels from Denver’s Open Data Portal (Existing Landuse 2018) and measured their area. A subset of these parcels near Union Station is mapped below.
Then, from all parcels, we sampled 100 sites of size greater than 1 acre, 25 sites of size 0.5-1 acre, and 400 sites less than 0.5 acres. One such sample appears below.
From this point, we followed the same framework as for other scenarios. The baseline NDVI of the sampled parcels and that of their corresponding 500m buffers is visualized below.
Finally, we propose replacing a portion of parking-lot surface with native plants. Denver’s area is about 9% parking lot, so even replacing a small amount of parking-lot surface could have a large impact on the total area of urban greening. We obtained spatial data on parking lots](https://www.denvergov.org/opendata/dataset/parking-lots-2016)) from the City of Denver Open Data Catalog.
| Total area (mi2), parking lots | Total area (mi2), Denver | proportion parking lots |
|---|---|---|
| 13.57 | 155 | 0.09 |
We similarly measured the NDVI on the parking lots as well as the NDVI in the census block groups within a 500-m buffer radius of any parking lot (which is almost the whole city). We present a small subset of these areas near Union Station:
We estimated the number of premature deaths averted under each scenario by following a recent meta-analysis that estimated that for every 0.1 unit increase in exposure to NDVI, the relative risk of premature death decreases about 4% (pooled risk ratio of 0.96; 95% confidence interval [CI]: 0.94, 0.97).3 We used this risk ratio to estimate the population attributable fraction9 corresponding to the proposed changes in NDVI from the baseline levels for each scenario to the alternative native-plants level.
Under the first two scenarios, the areal unit is the block group itself, and because we do not specify where, exactly, in the block group the greening intervention would occur, we define NDVI exposure as the average NDVI in the block group and compare that value with the alternative native-plants value. We use 5-year American Community Survey Data (2016-2020) to estimate the population in 5-year age groups in each census block group. For the other scenarios, we draw a 500-m buffer around the intervention area, as shown above. To estimate the population affected, we multiply the population density of the block group in that age group by the area of the block group covered by the intervention area.
To estimate the number of premature deaths prevented in each age group and area, we multiply the population-attributable fraction by the baseline mortality rate by the population size. We restrict analyses to adults aged 30 and above (subject to change) following the age range of many of the cohort studies reviewed3 and exclude northeastern census block groups near the airport (subject to change; perhaps not for all scenarios).
We estimate uncertainty in both the block-group-level population estimates and the relative risk estimate (95% CI: 0.94, 0.97) by resampling from normal distributions and taking the 2.5th and 97.5th percentiles over 1,000 replications. This code has additional detail.
Results of each scenario appear in the table below. The first table describes, for each scenario, the population affected, the cumulative area of the treatment itself (e.g., the area of the retention ponds or parking lots) and the cumulative area of the 500 m residential buffer surrounding the treatment area. Note for Scenarios 1 and 2, the treatment area and the residential-buffer area are the same (the area of the affected census block groups) because we do not say exactly where in the block groups the greening treatment would occur.
| Scenario | Sub-scenario | NDVI, native | Area, treatment (mi2) | Area, residential buffer (mi2) | Pop. affected, est. | Pop. affected, 95% LL | Pop. affected, 95% UL |
|---|---|---|---|---|---|---|---|
| all-bg | 100-pct | 0.34 | 1,161 | 1,161 | 152,144 | 147,306 | 156,154 |
| all-bg | 100-pct | 0.54 | 2,619 | 2,619 | 385,872 | 379,212 | 392,108 |
| all-bg | 30-pct | 0.34 | 1,161 | 1,161 | 152,144 | 147,604 | 155,841 |
| all-bg | 30-pct | 0.54 | 2,619 | 2,619 | 385,872 | 379,643 | 392,435 |
| all-bg | 20-pct | 0.34 | 1,161 | 1,161 | 152,144 | 147,777 | 155,912 |
| all-bg | 20-pct | 0.54 | 2,619 | 2,619 | 385,872 | 378,619 | 392,621 |
| riparian | 200-ft | 0.34 | 72 | 641 | 70,043 | 67,510 | 72,489 |
| riparian | 200-ft | 0.54 | 163 | 1,364 | 182,229 | 177,204 | 185,601 |
| riparian | 100-ft | 0.34 | 35 | 641 | 70,043 | 67,315 | 72,448 |
| riparian | 100-ft | 0.54 | 80 | 1,364 | 182,229 | 177,835 | 185,798 |
| riparian | 50-ft | 0.34 | 17 | 641 | 70,043 | 67,606 | 72,402 |
| riparian | 50-ft | 0.54 | 39 | 1,364 | 182,229 | 177,890 | 186,376 |
| ogi | ogi_proj | 0.34 | 11 | 261 | 45,386 | 43,003 | 47,623 |
| ogi | ogi_proj | 0.54 | 17 | 442 | 75,217 | 72,876 | 77,551 |
| ogi | parcel | 0.34 | 1 | 902 | 125,366 | 121,445 | 128,239 |
| ogi | parcel | 0.54 | 2 | 2,064 | 316,040 | 311,701 | 322,658 |
| prkng | 100-pct | 0.34 | 18 | 1,166 | 152,955 | 149,072 | 157,745 |
| prkng | 100-pct | 0.54 | 25 | 2,585 | 379,132 | 373,096 | 384,601 |
| prkng | 50-pct | 0.34 | 18 | 1,166 | 152,955 | 149,021 | 157,658 |
| prkng | 50-pct | 0.54 | 25 | 2,585 | 379,132 | 373,666 | 385,071 |
| prkng | 20-pct | 0.34 | 18 | 1,166 | 152,955 | 149,573 | 157,481 |
| prkng | 20-pct | 0.54 | 25 | 2,585 | 379,132 | 371,408 | 384,886 |
| all-bg-equity | cdphe-def-100-pct-20-pct | 0.34 | 400 | 400 | 36,066 | NA | NA |
| all-bg-equity | cdphe-def-100-pct-20-pct | 0.54 | 538 | 538 | 61,145 | NA | NA |
| all-bg-equity | cdphe-def-100-pct-30-pct | 0.34 | 536 | 536 | 56,956 | NA | NA |
| all-bg-equity | cdphe-def-100-pct-30-pct | 0.54 | 792 | 792 | 104,970 | NA | NA |
| all-bg-equity | cdphe-def-20-pct-20-pct | 0.34 | 400 | 400 | 36,066 | NA | NA |
| all-bg-equity | cdphe-def-20-pct-20-pct | 0.54 | 538 | 538 | 61,145 | NA | NA |
| all-bg-equity | cdphe-def-30-pct-30-pct | 0.34 | 536 | 536 | 56,956 | NA | NA |
| all-bg-equity | cdphe-def-30-pct-30-pct | 0.54 | 792 | 792 | 104,970 | NA | NA |
| all-bg-equity | denver-def-100-pct-20-pct | 0.34 | 404 | 404 | 34,069 | NA | NA |
| all-bg-equity | denver-def-100-pct-20-pct | 0.54 | 538 | 538 | 60,046 | NA | NA |
| all-bg-equity | denver-def-100-pct-30-pct | 0.34 | 494 | 494 | 48,559 | NA | NA |
| all-bg-equity | denver-def-100-pct-30-pct | 0.54 | 804 | 804 | 107,498 | NA | NA |
| all-bg-equity | denver-def-20-pct-20-pct | 0.34 | 404 | 404 | 34,069 | NA | NA |
| all-bg-equity | denver-def-20-pct-20-pct | 0.54 | 538 | 538 | 60,046 | NA | NA |
| all-bg-equity | denver-def-30-pct-30-pct | 0.34 | 494 | 494 | 48,559 | NA | NA |
| all-bg-equity | denver-def-30-pct-30-pct | 0.54 | 804 | 804 | 107,498 | NA | NA |
The second table shows, for each scenario, the estimated number of deaths prevented and the number of deaths prevented per 100,000 people affected.
| Scenario | Sub-scenario | NDVI, native | Deaths prevented, est. | Deaths prevented, 95% LL | Deaths prevented, 95% UL | Deaths prevented per 100k, est. | Deaths prevented per 100k, 95% LL | Deaths prevented per 100k, 95% UL |
|---|---|---|---|---|---|---|---|---|
| all-bg | 100-pct | 0.34 | 51 | 42 | 59 | 34 | 28 | 39 |
| all-bg | 100-pct | 0.54 | 304 | 288 | 328 | 79 | 75 | 85 |
| all-bg | 30-pct | 0.34 | 15 | 12 | 18 | 10 | 8 | 11 |
| all-bg | 30-pct | 0.54 | 88 | 82 | 95 | 23 | 21 | 25 |
| all-bg | 20-pct | 0.34 | 10 | 8 | 12 | 7 | 5 | 8 |
| all-bg | 20-pct | 0.54 | 59 | 54 | 62 | 15 | 14 | 16 |
| riparian | 200-ft | 0.34 | 1 | 1 | 2 | 2 | 1 | 2 |
| riparian | 200-ft | 0.54 | 12 | 10 | 13 | 6 | 5 | 7 |
| riparian | 100-ft | 0.34 | 0 | 0 | 0 | 0 | 0 | 1 |
| riparian | 100-ft | 0.54 | 4 | 4 | 5 | 2 | 2 | 3 |
| riparian | 50-ft | 0.34 | 0 | 0 | 0 | 0 | 0 | 0 |
| riparian | 50-ft | 0.54 | 2 | 1 | 2 | 1 | 1 | 1 |
| ogi | ogi_proj | 0.34 | 0 | 0 | 1 | 1 | 0 | 2 |
| ogi | ogi_proj | 0.54 | 2 | 1 | 3 | 2 | 1 | 3 |
| ogi | parcel | 0.34 | 0 | 0 | 0 | 0 | 0 | 0 |
| ogi | parcel | 0.54 | 1 | 1 | 1 | 0 | 0 | 0 |
| prkng | 100-pct | 0.34 | 15 | 13 | 17 | 10 | 9 | 11 |
| prkng | 100-pct | 0.54 | 48 | 43 | 51 | 13 | 11 | 13 |
| prkng | 50-pct | 0.34 | 7 | 7 | 8 | 5 | 4 | 5 |
| prkng | 50-pct | 0.54 | 24 | 22 | 26 | 6 | 6 | 7 |
| prkng | 20-pct | 0.34 | 3 | 3 | 3 | 2 | 2 | 2 |
| prkng | 20-pct | 0.54 | 9 | 9 | 10 | 2 | 2 | 3 |
| all-bg-equity | cdphe-def-100-pct-20-pct | 0.34 | 2 | NA | NA | 6 | NA | NA |
| all-bg-equity | cdphe-def-100-pct-20-pct | 0.54 | 11 | NA | NA | 18 | NA | NA |
| all-bg-equity | cdphe-def-100-pct-30-pct | 0.34 | 6 | NA | NA | 10 | NA | NA |
| all-bg-equity | cdphe-def-100-pct-30-pct | 0.54 | 29 | NA | NA | 27 | NA | NA |
| all-bg-equity | cdphe-def-20-pct-20-pct | 0.34 | 2 | NA | NA | 6 | NA | NA |
| all-bg-equity | cdphe-def-20-pct-20-pct | 0.54 | 11 | NA | NA | 18 | NA | NA |
| all-bg-equity | cdphe-def-30-pct-30-pct | 0.34 | 6 | NA | NA | 10 | NA | NA |
| all-bg-equity | cdphe-def-30-pct-30-pct | 0.54 | 29 | NA | NA | 27 | NA | NA |
| all-bg-equity | denver-def-100-pct-20-pct | 0.34 | 2 | NA | NA | 6 | NA | NA |
| all-bg-equity | denver-def-100-pct-20-pct | 0.54 | 9 | NA | NA | 15 | NA | NA |
| all-bg-equity | denver-def-100-pct-30-pct | 0.34 | 4 | NA | NA | 8 | NA | NA |
| all-bg-equity | denver-def-100-pct-30-pct | 0.54 | 25 | NA | NA | 23 | NA | NA |
| all-bg-equity | denver-def-20-pct-20-pct | 0.34 | 2 | NA | NA | 6 | NA | NA |
| all-bg-equity | denver-def-20-pct-20-pct | 0.54 | 9 | NA | NA | 15 | NA | NA |
| all-bg-equity | denver-def-30-pct-30-pct | 0.34 | 4 | NA | NA | 8 | NA | NA |
| all-bg-equity | denver-def-30-pct-30-pct | 0.54 | 25 | NA | NA | 23 | NA | NA |
As of April 27, 2022:
Add green streets HIA
Add economic-impact attributable to health
Copyright © 2022 Michael D. Garber